METHOD FOR MODELING SEMICONDUCTOR DEVICE PROCESS 



BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a method for 
modeling a semiconductor device process, particularly to a 
modeling method of impurity diffusion in a semiconductor, and 
a reverse snort channel effect of a threshold voltage of a 
MOS type field-effect transistor (MOSFET) . 

This application is a counterpart application of 
Japanese Application Serial Number 246268/2001, filed May 29, 
2001, the subject matter of which is incorporated herein by 
reference . 

Description of the Related Art 

FIG. 9 shows a mechanism of impurity diffusion in a 
case in which an excessive point defect exists in a 
semiconductor. An excessive point defect I is generated 
mainly in a high-dosage ion-implantation process. 
Particularly, the excessive point defect I generated in a 
source and drain (S/D) ion-implantation process of the MOSFET 
forms a pair with an impurity A by an immediately subsequent 
thermal treatment for allowing an ion to diffuse, and a pair 
AI of the impurity and point defect repeats separation and 
coupling, and diffuses. Thereafter, pairs AV and AI reach an 
interface between Si and SiOa, the point defect I of the pair 
disappears, and the impurity A having formed the pair piles 
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up. Additionally, V denotes a void, IV denotes a pair of the 
point defect and void, and AV denotes a pair of the impurity 
and void. 

FIG. 10 shows a relation between a distance of a 
5 channel direction Distance, and an impurity density Cone, in 
a case in which a gate length is set to 2.03 |a,m, 0.52 \ixa, and 
0.21 \m. In the drawing, a portion having a low impurity 
density corresponds to a channel, and a rise of the impurity 
density by impurity pileup is seen in portions corresponding 
IM to the drain and source on opposite sides of the channel. 

ip A pileup amount is largest in a gate end close to an 

iil 

-4; implantation position, and decreases toward the channel. 

hh 

"•-4 ^ reverse short channel effect of a threshold voltage of the 

MOSFET (a phenomenon in which the threshold voltage increases 
IS with a shorter gate length) is deemed to occur by the pileup 
=«' as a main factor. 

^ FIG. 11 shows a relation between a gate length Lg 

and a threshold voltage Vth when a substrate bias VB is OV, - 
3V, -5V. In a semiconductor device, when a device dimension 
20 is reduced so as to enhance a response speed and integration 
degree, the gate length is generally and accordingly 
shortened. 

On the other hand, the threshold voltage Vth is 
preferably constant regardless of the gate length Lg in a 
25 circuit design. However, the reverse short channel effect is 
seen. That is, when the gate length Lg is shortened (about 1 
lim or less in FIG. 11), the threshold voltage Vth vertically 
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fluctuates . 

A conventional model for describing impurity 
diffusion in the semiconductor is roughly classified in three. 
That is, a first model is a method including: setting an 
5 analysis object region; solving one diffusion equation with 
respect to each impurity; and calculating an impurity 
distribution in the semiconductor (hereinafter referred to as 
a Fair model) . A second model is a method including: setting 
the analysis object region; assuming that the point defect 

]|,0 and impurity form a pair and diffuse, and solving respective 
J« diffusion equations with respect to the point defect, and the 
pair of impurity and point defect; and calculating the 

V impurity distribution in the semiconductor (hereinafter 

referred to as a pair diffusion model) . A third model is a 

I'B simple model for setting an analysis region, and simply 
taking the effect in a frame of the Fair model without 

,5^: solving the diffusion equation associated with the point 
defect. 

The following equations are equations for use in 
20 using the Fair model and pair diffusion model to obtain the 
impurity distribution in the semiconductor. 
Fair model 

dCa/dt = -VJa (Equation 12) 
Pair diffusion model 
25 3Ci/dt+i:aCAi/3t = -VJi-I:VJai-Ri,v+Riv (Equation 13) 

aCv/at+EaCAv/3t = -VJv-EVJav-Ri,v+Riv (Equation 14) 

aCAtotal/at = V{DAiCAtotal(Ci/Ci*)V(log(pCAtotalCi/Ci*ni)} 
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+V{ (Cv/Cv* ) V ( log (pCAtotaiCv/Cv*ni) } 

(Equation 15) 



10 

m 
m 
m 

m 



density (cm"^) Cfttotai = ICa+ICai+ICa^ 
flux (cm"^/s) 

re-coupling speed (cm"-^/s) 
carrier density (cm"^) 
equilibrium density (cm"'^) 



ni: intrinsic carrier density (cm" ) 
D: diffusion constant 
A: all dopant species 
I : interstitial 
V: vacancy 

AI: pair of dopant species and interstitial 
IV: pair of interstitial and vacancy 
AV: pair of dopant species and vacancy 
As shown above, for the first Fair model, since the 
number of equations to be solved is small, a calculation time 
is short, and the model can advantageously easily be handled. 
On the other hand, an influence of the point defect with 
respect to the impurity diffusion is not handled, and 
therefore there is a disadvantage that the impurity pileup in 
the interface of Si and Si02 cannot be reproduced. 

In the second pair diffusion model, since the 
influence of the point defect with respect to the impurity 
diffusion is strictly handled, simulation can advantageously 
be performed with a high precision. On the other hand, when 



the number of impurities increases, the number of equations 
to be solved also increases, and the calculation time 
disadvantageously increases. 

The third simple model has an advantage that 
tradeoff of the calculation time and simulation precision can 
well be absorbed, but a degree of reproduction of an 
application range of the model and process dependence is a 
key. 

Examples of the simple model include a simple model 
disclosed, for example, in Japanese Patent Application Laid- 
Open No. 084716/2000. FIG- 12 is an explanatory view of the 
simple model, FIG. 13 is a diagram showing a lateral 
distribution of the impurity pileup in an end of a gate 
electrode in the semiconductor substrate, and FIG. 14 is a 
flowchart showing the simple model. 

In FIG. 12, C(xi,yj) denotes the impurity density of 
each cell, AC(xi,yj) denotes the impurity density which 
decreases in each cell, and AS(xi,yj) denotes an area of each 
cell. An amount of the impurity included in the cell (xl,yl) 
at the interface is C (xl, yl) AS (xl, yl) before the diffusion, 
and C(xl,yl)AS(xl,yl)+AC(xi,yj)AS(xi,yj) after the diffusion. 
In FIG. 13, Cint(x) denotes the impurity pileup amount in the 
interface, and Cpile denotes the impurity pileup amount in a 
region in which a gate 10 is not disposed. FIG. 13 shows 
that the pileup amount of the impurity decreases according to 
Equation 16 in a portion closer to a middle portion of the 
gate 10. Moreover, the impurity pileup amount is constant at 



Cpile in the region where the gate 10 is not disposed. 

Cint(x) = Cpile exp (-x/A.int ) {Equation 16) 

In FIG. 14, t denotes time. Particularly, tiED 
denotes a transient enhanced diffusion (TED) duration. 

As shown in FIGS. 12, 13 and 14, first an Si layer 
is divided into a plurality of cells. For the impurity 
pileup in the interface of Si and Si02 which cannot be 
reproduced by the Fair model, instead of solving the 
diffusion equation associated with the point defect, a part 
of an impurity amount of a substrate region is moved to the 
interface of Si and Si02 and constitutes the pileup. An 
impurity movement mass is used as a function of a distance of 
a noted interface position to the cell of the substrate 
region in a method for use. 

When the diffusion model is used to model the 
reverse short channel effect of the threshold voltage, for 
example, the following methods are used: 

(1) a method of using the pair diffusion model to 
calculate the impurity distribution, and calculating an 
electric property in an unchanged state; and 

(2) a method of using the simple model to calculate 
the impurity distribution, and calculating the electric 
property in the unchanged state . 

However, the aforementioned methods have some 
problems . 

1. In the pair diffusion model, the calculation 
time increases. Therefore, when simulations such as 
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sensitivity analysis among the process;- device, and circuit, 
process optimization, process dispersion analysis, and 
calibration of a model parameter need to be executed a 
plurality of times, it is difficult to use the pair diffusion 
5 model. 

2. In the conventional simple model, since the 
calculation of a mass of impurity moving to the interface of 
Si and Si02 is represented by the function only of the 
interface position and the distance to the cell of the 

10 substrate region, dependence of the S/D process on an 
impurity density re-distribution cannot be reproduced. 

"|| 3. Moreover, in the conventional simple model, a 

N' shape of the impurity pileup in the interface is determined 
beforehand. Therefore, there is a problem that the 

%^ dependence of the process on the shape of the pileup cannot 

y be predicted. 

m 

SUMMARY OF THE INVENTION 
An object of the present invention is to provide a 
20 method for modeling a semiconductor device process in which: 
S/D process dependence of an impurity pileup in an 
interface between Si and Si02 can be reproduced in a frame of 
a Fair model, 

a shape of the impurity pileup can be predicted, 
25 an impurity profile can be used to reproduce a 

reverse short channel effect of a threshold voltage, and 
a high-speed calculation is possible. 
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To achieve the object, according to the present 
invention, there is provided a method comprising: a first 
step of setting data of an Si02 layer; a second step of 
setting data of an Si layer formed in contact with the Si02 
layer; a third step of setting a plurality of cells in the Si 
layer, and setting an amount of an impurity included in each 
cell; a fourth step of setting an amount per unit time by 
which the impurity included in each cell moves to another 
cell; a fifth step of setting data by which a cell in the 
vicinity of an interface of the SiOa layer and the Si layer 
is set as an impurity pileup portion; a sixth step of setting 
data of a position of a source or a drain in the Si layer; 
and a seventh step of calculating the amount of the impurity 
included in each cell for each unit time after each of the 
first step, the second step, the third step, the fourth step, 
the fifth step, and the sixth step ends, wherein a mass of 
impurity moving to the pileup portion from each cell is 
determined as a function of a distance to the impuriry pileup 
portion from each cell (hereinafter referred to as distance 
rl), and a distance to the source or the drain from each cell 
(hereinafter referred to as distance r2) . 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 is a block diagram showing a computer for 

simulation and data to be inputted into the computer. 

FIG. 2 is an explanatory view of a method for 

modeling a semiconductor device process according to a first 
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embodiment . 

FIG. 3 is an explanatory view of the method for 
modeling the semiconductor device process according to a 
second embodiment . 
5 FIG. 4 is an explanatory view of the method for 

modeling the semiconductor device process according to a 
third embodiment - 

FIG. 5 is an explanatory view of the method for 
modeling the semiconductor device process according to a 
:|0 fourth embodiment . 

FIG. 6 is an explanatory view of the method for 
modeling the semiconductor device process according to a 
N= fifth embodiment. 

rt*^ FIG. 7 is a flowchart showing a process of the 

i;^- method for modeling the semiconductor device process 
y according to the fifth embodiment. 

'^l FIG. 8 is an explanatory view of the method for 

iy 

modeling the semiconductor device process according to a 

sixth embodiment. 
20 FIG. 9 is an explanatory view of a mechanism of 

impurity diffusion. 

FIG. 10 shows a characteristic curve of a relation 

between a channel longitudinal direction distance and an 

impurity density. 
25 FIG. 11 shows a characteristic curve of a relation 

between a gate length and a threshold voltage. 

FIG. 12 is an explanatory view of a conventional 
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s-iraple model . 

FIG. 13 is a diagram showing a distribution of a 
pileup amount in an Si/SiOa interface. 

FIG. 14 is a flowchart showing a process of the 
5 method for modeling the semiconductor device process in the 
conventional simple model . 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Preferred embodiments of a method for modeling a 
JO semiconductor device process according to the present 
invention will be described hereinafter in detail with 

hi 

yj reference to the accompanying drawings. Additionally, in the 
bl' present specification and drawings, constituting elements 
111 having substantially the same functions/constitutions are 
1|| denoted with the same reference numerals and redundant 
description thereof is omitted. 
(First Embodiment) 

In a method for modeling the semiconductor device 
process according to a first embodiment, for a conventional 
20 pimple model which cannot reproduce S/D process dependence, 
instead of solving a diffusion equation associated with a 
point defect, a part of an impurity amount of a substrate 
region is moved to an interface of Si and Si02 to constitute 
an impurity pileup, and an impurity movement mass is 
25 determined as a function of a distance rl from a pileup 

position and a distance r2 to an S/D profile position of a 
closest interface. 
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In the conventional simple model, for example, when 

an S/D implantation dosage is changed, dependence of an S/D 

process on the impurity pileup in the interface cannot be 

reproduced. This is because information on the S/D process 

5 is not reflected with respect to calculation of the mass of 

impurity moving to the interface of Si and Si02 from each 

cell in an Si layer. 

FIG. 1 is a block diagram showing a computer for 

simulation and data to be inputted into the computer 

10 according to the present invention, and FIG. 2 is an 

explanatory view of the method for modeling the semiconductor 

device process according to the first embodiment. As shown 

H= in FIG. 1, a gate 112 is formed on an upper Si02 layer 110 of 

f|'| a semiconductor element 100 via an interface of Si and SiOa 

1^, 130, and a drain 126 is formed on a lower Si layer 120. An 

y ' ^^'^ profile position 132 changes with conditions of 
M= 

implantation, and the like. 
^ It is now assumed that the Si layer 120 is divided 

into a plurality of cells, and an impurity pileup amount 
20 Cint(x) in a cell 122 of the interface is calculated based on 
the following Equation 1 . 

Cint(x) = EAC(xi,yj) (Equation 1) 

It is assumed that an original impurity amount 
C(xi,yi) moves from an arbitrary cell 124 having an area 
25 AS(xi,yi) and an impurity amount AC(xi,yi) piles up in the 
cell 122. As shown in FIG. 1, first, a step of inputting 
data of Si02 and Si layers into the computer is performed. 
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That is, data concerning the S/D process dependence and Si 
layer density dependence, and the like are inputted. 

An impurity movement mass AC(xi,yi) is determined as 
the function of the distance rl from the pileup position and 
5 the distance r2 from the S/D profile position 132 of the 

closest interface, so that the S/D process dependence of the 
mass of impurity moving to the interface 130 can be 
reproduced based on the following Equations 2, 3, and 4. 

AC(xi,yj) = AXPbxPixC(xi,yj)xAS (xi,yj) (Equation 2) 
ig Pb = exp(-rl/^l) (Equation 3) 

|ij Pi = exp(-r2/;i2) (Equation 4) 

Parameters: A, a, B, Xl, X2 
C(x,y): impurity density in (x,y) 
\P- AS(x,y): cell area in (x,y) 

]|S Here, a parameter A relates to a magnitude of a 

U| reverse short channel effect, a and A-1 are parameters 

w 

concerning the S/D process dependence, and p and X2 are 
* parameters concerning Si layer density dependence. 

As an example, calculation can be performed assuming 
20 that a size of the cell AS is 10~^|jm^, original impurity 
amount C is lO^^cm"^, and impurity movement mass AC is 
2xl0^^cm"^ (20% of the original impurity amount C) . Moreover, 
it can be assumed that the parameter A, S/D process 
dependence a and Xl, and parameters p and X2 concerning the 
25 Si density dependence are 3.7x10"^, 1.0, 2.0 ^m, 1.0, 0.5 )jm, 
respectively. 

For example, when the S/D implantation dosage 



increases, the S/D junction position 132 of the interface 
shifts to a channel side, therefore the distance r2 is 
shortened with respect to the cell in the vicinity of the 
channel, the mass of impurity moving to the interface 130 
increases, and as a result the S/D process dependence of the 
interface on the impurity pileup can be reflected. 

According to the first embodiment, the S/D profile 
dependence onto the impurity pileup in the interface can 
easily be taken, and it is further possible to predict a 
pileup shape with respect to a process condition. 
(Second Embodiment) 

In the method for modeling the semiconductor device 
process according to a second embodiment, an S/D profile 
reference position of the distance r2 is determined as a 
distance from a constant density line 127 of an S/D profile 
in the first embodiment. 

FIG. 3 is an explanatory view of the method for 
modeling the semiconductor device process according to a 
second embodiment. A size and shape of the impurity pileup 
in the interface 130 depend on an excessive point defect 
density generated by an S/D implantation damage, and also 
influence an initial point defect distribution. 

The initial point defect distribution is regarded as 
an S/D shape, and the impurity movement mass is determined as 
a. function of a distance rl from the pileup position and a 
distance r2 from a closest constant density line 127 of the 
S/D profile, so that the S/D process dependence of the mass 
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pf impurity moving to the interface can be reproduced. 
Equations for use conform to those of the first embodiment. 

According to the second embodiment, there is an 
advantage that an effect of the S/D shape changing in 
5 accordance with an energy and dosage of S/D can be reflected 
in impurity .pileup calculation in the interface 130. 
(Third Embodiment) 

In the method for modeling the semiconductor device 
process according to a third embodiment, a solid angle is 
\2 considered in calculating the mass of impurity moving to the 
interface 130 from the arbitrary cell 122 in the Si layer 
such as the substrate in the first and second embodiments. 
FIG. 4 is an explanatory view of the method for 
=f modeling the semiconductor device process according to the 
15 third embodiment. In the calculation of the mass of impurity 
^, moving to the interface 130 from the arbitrary cell 122 in 
<n the Si layer, even when the distance rl to the cell 124 of 

the interface from the arbitrary cell 122 in the Si layer is 
the same, but when a sectional area obtained by connecting 
20 the cell 124 of the interface to the arbitrary cell 122 in 

the Si layer differs, then the mass of impurity moving to the 
interface also differs. Therefore, a solid angle A9 is 
considered in the calculation of the mass of impurity moving 
to the interface from the arbitrary cell 122 in the Si layer. 
25 The equations for use conform to those of the first 

embodiment, but Equation 6 is obtained by inserting a term of 
parameter Pr concerning the solid angle into Equation 2. 



- 15 - 



Equation 6 is used instead of Equation 2. 

AC(xi,yj)=AXPbxPixPrxC(xi,yj)xAS(xi,yj) (Equation 6) 
Arbitrary cell in the Si layer 
The following method is an example of a calculation 
5 method of the solid angle A9. First, a side (a) of the 

interface cell 124 in the interface, side (b) connecting one 
of two ends of the side (a) to a representative point of the 
cell 122, and side (c) connecting the other end of the side 
(a) to the representative point of the cell 122 are 
10 calculated. 

5|l Subsequently, the sides (a), (b) , and (c) are used 

y-| to calculate the solid angle A6 by the following equations. 

'§ 

l# Ae = ATAN(r/ (p-a) 

tyi p = 0.5* (a+b+c) 

l| r = rl"(l/2) 

y - rl = (p-a) * (p-b) * (p-c) /p 



(*quoted from p. 174 of ^'Basic Mathematics Handbook" 
authored by Miyamoto, et al . , published by Morikita Shuppan, 
1990) 

This A6 is assigned to Equation 7, and Equations 1, 
3, 4, and 6 are used to calculate the impurity pileup amount 
Cint (x) . 

Pr = Ae/27C (Equation 7) 

According to the third embodiment, in the 
calculation of the mass of impurity moving to the interface 
from the arbitrary cell in the Si layer, since the solid 
angle is considered, position dependence between two cells 
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can advantageously be reflected in forming the impurity 
pileup in the interface. 
(Fourth Embodiment) 

In the method for modeling the semiconductor device 
5 process according to a fourth embodiment, the calculation of 
the impurity movement to the interface from the arbitrary 
cell in the Si layer is performed also using generation and 
disappearance terms of the diffusion equation in the first, 
second, and third embodiments. The equations for use conform 

10 to those of the first embodiment. 

FIG. 5 is an explanatory view of the method for 

'oi modeling the semiconductor device process according to the 

fourth embodiment. In the conventional simple model, for the 
impurity movement to the interface from the arbitrary cell 
122 in the Si layer, it is presumed that the impurity 
momentarily moves to the interface out of the Si layer. This 
is because the point defect diffuses very fast. Additionally, 
in fact, during movement, some of impurities remain in the 
substrate, and this influence becomes clear in a lower- 

20 temperature heat treatment. 

Therefore, some of the movement mass of impurities 
are assigned to the impurity generation and disappearance 
terms of the diffusion equations such as Equations 7 and 8, 
and a method of direct addition to the cell of the interface 

25 is also used, so that the impurity on the cell being moved is 
indirectly influenced. This is a characteristic of the 
method for modeling the semiconductor device process 
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according to the fourth embodiment- Here, since a total mass 

of impurities has to be stored in an analysis region, 

generation and disappearance amounts of impurities have to be 

set to be equal to each other. 
dC 

— +divJ = 0 (Equation 7) 

dC 

— +divJ = -Ac* or +AC* (Equation 8) 
at 

Generation: +AC* 
Disappearance: -AC* [cm"^/s] 

Here, 

J = -D grad C (Equation 9) 

C: impurity density 
D: diffusion constant 
t: time 

AC*: generation or disappearance amount of C 

As shown in FIG. 5, some of the mass of impurities 
are allowed to disappear in the arbitrary cell 122 in the Si 
layer based on Equation 7 (-AC*), and generated on the cell 
124 of the interface (+AC*) . This indirectly influences the 
respective cells through a flow of impurities FLAX (J) shown 
by arrow marks in accordance with Equation 8 . 

As an example, calculation can be carried out 
assuming that -AC* is 2 X 10^^ cm"^ and +AC* is 2 X 10^^ cm~\ 

According to the fourth embodiment, an effect that 
some of the mass of impurities moving to the interface from 
the respective cells in the Si layer remain during movement 
can indirectly be taken, and an application range can 
advantageously be extended also to the low- temperature heat 
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treatment in which the effect becomes clear. 
(Fifth Embodiment) 

In the method for modeling the semiconductor device 
process according to a fifth embodiment, in the first and 
second embodiments, an exceptional processing is performed in 
a case in which a plurality of S/D regions and isolated 
regions formed of the same type of impurities exist. In the 
method, a mask is used to segment a region formed as an 
object by implantation, and a region formed beforehand. 

FIG. 6 is an explanatory view of the method for 
modeling the semiconductor device process according to the 
fifth embodiment, and FIG. 7 is a flowchart showing the 
method for modeling the semiconductor device process 
according to the fifth embodiment. 

As shown in FIG. 6, in a semiconductor element 500, 
except an implantation region contributing to the impurity 
pileup in the interface, an isolated region 532 formed of the 
same type of impurities beforehand exists in the distance r2 
to an S/D region position from each cell in the Si layer 
treated in calculating the movement mass of impurity moving 
to the interface. In this case, a mask 505 is used, and the 
isolated region is excluded from a search object. 
Additionally, the S/D region mentioned herein refers to a 
region in which an excessive point defect is generated by 
implantation, and does not necessarily refer to the S/D 
region. 

In the example of FIG. 6, first an implantation 



region formed of the same type of isolated impurities is 
searched in an analysis region excluding a range defined by 
the mask. For example, a cell 522 refers to a closest region 
530, and a cell 524 refers to a closest region 528. Since 
the closest region 532 of a cell 526 is in the range defined 
by the mask 505 and is other than the search object, the 
region 528 is referred to. Additionally, reference numeral 
533 denotes an element isolation region, and the cells in the 
Si layer existing on the right side of the element isolation 
region 533 are excluded from the search object in the fifth 
embodiment . 

As shown in FIG. 7, based on the aforementioned 
method, first the isolated implantation region is searched in 
the analysis region outside the range defined by the mask 
(step 601), and the mass of impurities moving to the 
interface from respective cells such as the cells 522, 524, 
and 526 is calculated and moved (step 602) . 

Subseguently, the diffusion equation is solved by 
the Fair model (step 603) . When a heat treatment time t does 
not exceed a time tTED (transient enhanced diffusion duration) 
(step 604), the flow returns to the step 601. When the heat 
treatment time exceeds the time tTEo^ the flow advances to a 
normal diffusion process. 

According to the fifth embodiment, when a mask 
designation function already mounted usually on software is 
used, the region formed as the object by the implantation can 
advantageously and easily be segmented from the other regions 
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at a small development cost. 
(Sixth Embodiment) 

In the method for modeling the semiconductor device 
process according to a sixth embodiment, an impurity profile 
5 constituted in the first to fifth embodiments is used to 
model the reverse short channel effect of the threshold 
voltage. The equations for use conform to those of the first 
embodiment . 

In general, parameters corresponding to the 
W respective processes are introduced in order to reproduce the 
process dependence of the reverse short channel effect of the 

\n 

\0 threshold voltage. In the sixth embodiment, Xl and 12 

associated with the distance rl to the interface cell from 
the arbitrary cell in the Si layer and distance r2 to the S/D 
region during calculation of the mass of impurity moving to 
the interface are used as a function of a diffusion length of 
the point defect to represent the following equations. 
A-1 = ax(DiXtTED) (Equation 10) 
X2 = j3x(DiXtTED) (Equation 11) 

Here, a and P are parameters, and a is a parameter 
representing the S/D process dependence, and absorbs the 
process dependence in the S/D formation. The parameter P 
represents the Si layer density dependence, and absorbs the 
process dependence such as' a dependence of a threshold 
voltage control implantation condition. Moreover, a 
parameter A for adjusting the magnitude of the reverse short 
channel effect of the threshold voltage is introduced. 
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In the sixth embodiment, since the parameter 
corresponding to the process condition is introduced in this 
manner, the process condition dependence of the reverse short 
channel effect of the threshold voltage can be reproduced, 
^nd experimental measurement is easily combined at a small 
calculation cost. 

FIG. 8 shows a result of simulation of gate length 
dependence of the threshold voltage by the calculation method 
together with experimental measurement values. Here, 
characteristic curves concerning three process condition Nos. 
1, 2, 3 are shown, and in any case a result remarkably 
approximate to the experimental measurement value is obtained. 

The preferred embodiments of the method for modeling 
the semiconductor device process and program according to the 
present invention have been described above with reference to 
the accompanying drawings, but the present invention is not 
limited to these embodiments. A person skilled in the art 
would apparently develop various modifications in a category 
of technical thoughts described in the appended claims, and 
such modifications are also naturally deemed to fall within a 
technical scope of the present invention. 

For example, in the method for modeling the 
semiconductor device process according to the present 
invention, manufacturing may be carried out by an SOI process. 
In this case, an Si layer sandwiched between a gate oxide 
film and a BOX layer, and an Si layer as the substrate exists - 

Among these Si layers, the impurity pileup by the 
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heat treatment is not generated in the Si layer as the 
substrate. This is because the impurity pileup is determined 
by a positional relation between an atmospheric air (gas and 
air) and respective layers. Therefore, in this case, when 
5 the parameter is selectively adjusted with respect to two 

interfaces between the gate oxide film and BOX layer, and the 
sandwiched Si layer, the method of the present invention can 
be applied. 

Moreover, a simulation part from when the process 
;3:0 condition is inputted until the electric property is 
f*II outputted can also be replaced with an apparatus. 

m 

■■C As described above, according to the present 

■y invention, there can be provided a method for modeling the 

m 

55 semiconductor device process in which the S/D process 

dependence of the impurity pileup in the interface of Si and 
l^h Si02 and the shape of the impurity pileup can be reproduced 
j^'i in the frame of the Fair model, the impurity profile can be 
used to reproduce the reverse short channel effect of the 
threshold voltage, and high-speed calculation is enabled. 



